“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1968 


A study of the Hermitian numerical method 
applied to the single degree of freedom, 
damped oscillator. 


Van Sickle, Garth Allan 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/13163 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 


1968 
VAN SICKLE, G. 





A STUDY (OF THE HERMITIAN NUMERICAL. 
METHOD > APPL [to TO THE SINGLE DEGREE OF FREEDOM, 
DAMPED OSCILLATOR | 
. GARTH ALLAN VAN SICKLE 














A STUDY OF THE HERMITIAN NUMERICAL METHOD APPLIED 
TO THE SINGLE DEGREE OF FREEDOM, DAMPED OSCILLATOR 


by 


Garth Allan Van Sickle 
Ensign, er gee Navy 
B.S., U.S. Naval Academy, 1967 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 
from the 
NAVAL POSTGRADUATE SCHOOL 
March 1968 


ABSTRACT 


Approximation of the differential equation of the single 
degree of freedom, damped oscillator by the first order central finite 
difference equations and two Hermitian finite difference equations is 
investigated. An error analysis is made between the solution to the 
differential equation and the solutions to the finite difference 
approximations for various values of damping. The results of the error 
analysis indicate that the Hermitian approximations have less error 
than the first order difference approximation for the same size of 
increment. Furthermore, the employment of the Hermitian method does 
not materially increase the execution time on the digital computer over 
that required by the first order difference equations. Thus the 
Hermitian equations are superior to the first order finite difference 


equations for the damped oscillator problem. 
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CHAPTER I 


INTRODUCTION 


In the quest for accurate numerical solutions to differential 
equations, considerable effort has been devoted to the determination 
of the errors associated with various finite difference schemes. This 
work has gained importance with the recent widespread use of the 
digital computer. Due to the high value placed on computer time and 
to the limitations on core storage space within the computer, the 
accuracy and efficiency of the numerical method employed are of the 
utmost importance. 

In general, many engineering problems can be classified as 
either boundary value problems or initial value Selene For the 
numerical solution to boundary value problems, the most common method 
used is a form of the central finite difference operators. Several 
forms of these operators can be found in Table III of the Appendix to 
reference 1. Some of these are frequently used to solve engineering 
problems while others seem either to be impractical or to appear 
seldom in current literature. The Hermitian equations, which are 
included in this table, fall into the last category. For initial 
value problems, numerical methods of solution such as Adams, Runge- 
Kutta, and other forms of the predictor-corrector methods are usually 
employed. The Hermitian method is also applicable to initial value 
problems but is seldom used. 

The purpose of this paper is to evaluate the merits of the 
Hermitian method when applied to a second order, ordinary, linear 
differential equation. This evaluation will be made by comparing the 
solution to the exact differential equation with the solution to two 
Hermitian finite difference approximations and with the solution to 
the first order central finite difference equation. In this manner, 
the difference in accuracy between these finite difference methods can 


be determined, and the merits of the Hermitian method established. 
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The equation chosen for this analysis is that of the single 
degree of freedom, damped oscillator. This equation was chosen | 
because it contains both the first and second derivatives and because 


+3 This appears to be the first | 


its exact solution is well known. 

application of the Hermitian equations to differential equations con- 

taining more than one degree of derivatives. This problem is generally 

considered an initial value problem. However, the results of this 

analysis are also applicable to boundary value problems since a boundary 

value can be solved using methods of solution normally employed in 

initial value problems. Also, the oscillator may become a boundary 

value problem if the displacement or velocity is specified at a time 

other than zero. Thus, the choice of this example allows the conclu- 

sions to be applied to both initial value and boundary value problems. | 
The author wishes to acknowledge a special thanks to 

Professor Robert E. Ball of the Aeronautics Department of the Naval | 

Postgraduate School. The guidance and consultation of Dr. Ball, as 


thesis advisor, were instrumental in the development of this paper. 
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CHAPTER II 


DISCUSSION OF THE CENTRAL FINITE 
DIFFERENCE METHOD AND ITS ASSOCIATED ERRORS 


There are three types of error associated with the finite 
difference solution to engineering puewlene One type is called 
truncation or discretization error. The conversion of a differential 
equation to an algebraic equation is, in essence, the approximation 
of a continuous system by a discrete mode. The error involved in 
this conversion is called the truncation error. A graphic example of 
how this error is introduced can be seen by examining the derivation 
of the first order central finite difference equations for the first 
and second derivatives presented in Appendix A. 

In the derivation of the finite difference expression for the 
first derivative, only the first three terms of the Taylor series are 
used. The first four terms are used in the expression for the second 
derivative. The truncation error can be approximated by the first 
unused term in the Taylor series. Therefore, a decrease in truncation 
error can be accomplished by using more terms of the Taylor series in 
the finite difference equation. 

Another type of error is called round-off error. This error 
is generated within the computer and is a function of the amount of 
calculation required and the condition of the numbers involved. Each 
arithmetic computation, such as addition, multiplication, and expo- 
nentiation, has a round-off error associated with it. The computer 
is limited to carrying a specified number of digits. Any amount over 
this number must be dropped for each calculation. To decrease round- 
off error, it is essential to minimize the number of arithmetic calcu- 
lations within the program and to avoid subtraction of large numbers, etc. 

The third type of error is called inherent error. This error, 
unlike truncation error and round-off error, is not caused by the use 
of numerical methods.” Inherent error is introduced into the program 
within the data. If the data are not exact, which is usually the case, 


the solution cannot be expected to be exact. 
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The object of any numerical analysis is to obtain as accurate 
an answer as possible without significantly increasing the computation 
time. There are a number of higher order finite difference methods 
available that have less truncation error, Hence, they prowide more 
accurate solutions than the first order central finite difference 
method, Among these is the Hermitian method. These methods vary by 
the number of stations used, the location of the stations, ete. 

When central finite difference methods are applied to an 
engineering problem, the problem is reduced to that of solving a 
matrix. This matrix will usually be a band matrix centered about the 
principal diagonal. The width of the band is dependent upon the number 
of stations included in the finite difference equation. Using the first 
order difference equations listed in Appendix A, the band width would 
be three, and the matrix is called a tridiagonal matrix. If higher 
order difference equations that require more stations are used, the 
band width widens. By widening the band of the matrix, more computa- 
tions may be required than for the tridiagonal matrix. This increases 
the round-off error and the computation time. Furthermore, the tri- 
diagonal matrix will require less core storage than the wider band 
matrix for the same element size. 

Another approach used to decrease the magnitude of the trun- 
cation error involves decreasing the interval between the stations. 
This is equivalent to increasing the degree of the matrix and, there- 
fore, introduces greater round-off error. Other effects of this pro- 
cedure are the increase in computer storage space required and the 
increase in computer running time. 

The Hermitian method considered in this thesis results in a 
tridiagonal matrix but also decreases the truncation error. Conse- 
quently, the gain in accuracy due to the decrease in truncation error 
is not offset by an increase in round-off error. Furthermore, this 
method can be used to advantage on those problems where core storage 


requirements are important. 
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CHAPTER III 


APPLICATION OF THE HERMITIAN METHOD 
TO THE DIFFERENTIAL EQUATION 


The first order central finite difference equations for the 


first and second derivative, as derived in Appendix A, are 


y= sly, , = Je, | aes) 








— ae 
A? Sys 29; 4)-1 | (3.1b) 
The Hermitian equations for the first and second derivatives, as 


derived in Appendix B, are 


[hes Yi SY 


ui 
3 \vt ay, + a] 





(3.2a) 


44 
[y.2y+¥)- 4 Aly Zr Hoy +Y 4", (3.2b) 
where (') denotes the oa, ci respect to time. 

As shown in Appendix B, the Hermitian equations have less trun- 
cation error than the first order equations. Note, that in order to 
get this reduction in trucation error it was not necessary to use 
more stations. 

Consider the single degree of freedom, damped oscillator shown 


in Figure 1. The differential equation of the oscillator is 459 


mys cys ky —0 (3.3) 
where m is the mass of the system, c is the viscous damping coeffi- 
cient, k is the spring constant, and y is the displacement. Sub- 
stituting the conventional first order central difference equations, 
Eqs. 3.la and 3.1b, for the first and second derivatives into Eq. 3.3 


leads to 


Bycey +Y),] y at Yar ¥,.,] j KY) eis Se. 


where /\ has been replaced with{\ ¢ to denote the time interval. 
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DAMPED OSCILLATOR 
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Since the Hermitian equations are not of the same form as 
the conventional equations, a different method of approximation to 
Eq. 3.3 is necessary. Examination of the right side of Eqs. 3.2a 
and 3.2b reveals that Eq. 3.2a contains the first derivative at the 
three stations j+l, j, and j-l. Likewise, the right side of Eq. 3.2b 
contains the second derivative at the same three stations. By adding 
a third equation expressing y at j+l, j, and j-1 and modifying one of 
the Hermitian equations, it is possible to solve the differential 
equations at these three stations simultaneously. Equations 3.5, 


3.6, and 3,7 show this procedure. 




















mn maak ae or 
eS =| | le) }= 12 (yj. #10 y: Fy | 
(345) 
a Je Te Waeayiegs, | 
4 At i445 — Ba See Ae 
(3.6) 
eae + — kk 
j2 Yi Oy To), |= |2 Y.,t10y, ae | 
(3.7) 
Adding Eqs. 3.5, 3.6, and 3.7 gives 
mM _ Cc a a 
Are (Ys oe } + Tar (Ya, ae | 
he _ (OM SaGas 7, | ON 
ape = yy, toy ty, | = jo a Ope pp gwetty 





(3.8) 


Li 


Note that the right side of Eq. 3.8 will equal zero if 1 y'; is 
added to both sides. The conventional approximation for yj gives 


C 7 


a Sy am jo. 4 








(3.9) 
If Eq. 3.9 is added to Eq. 3.8, the result is 
mM Cc 
At2 4, 2Y, ed | . 2 At be so 
+ (ye) Hoy ty } = O 
a ela id 
(3.10) 


This is one form of the Hermitian approximation to Eq. 3.3. 
To get the second form of the approximation, it is necessary 


to write Eq. 3.7 as 


[Yat 4Y Y)= Ely. +49 "Ye | 





(3.11) 
When Eqs. 3.5, 3.6, and 3.11 are added, then 
ar tYyeulat sere aor | 
Arydt Uv “Ua 2At\U+! J—! 
k — 10M 4 4ce 44k 
+ Ws *4Y, “y= PS J" 6% 
ie 


Note that the right side of Eq. 3.12 will equal zero if my 5 is 


subtracted from both sides. From the conventional approximation for 


ve 


le 


m 


oe 





Yi 2y ty, J=myy ean 
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If Eq. 3.13 is subtracted from Eq. 3.12, the result is 


mM 


= ge - 
ar lye &Y +y,.) ey sat Yn U- 
Bie = 
6 [Yuh 4y, +y, | : (3.14) 


This is the second form of the Hermitian approximation to Eq. 3.3. 
In this thesis, Eq. 3.10 will be called MOD 2, and Eq. 3.14 will be 
called MOD 1. 

By following the derivation of these two Hermitian approxi- 








mations, it should be clear that a first order central finite 
difference equation for either y', or yey was used in each deriva- 
tion. For example, the first order central finite difference equa- 
tion for y" 5? Eq. 3.9, was introduced in the derivation of Eq. 3.10, 
the MOD 2 approximation. This is equivalent to modifying Eq. 3.6 

to be 


Sly, —-y l= Sy, HOoyrty. 

DAU Ut J 12 ae Yd as) 
Appendix C shows that the truncation error in Eq. 3.15 is the same 
as the truncation error in the first order central finite difference 
equation for the first derivative, Eq. 3.9. Thus, when the Hermitian 
equations are applied to the differential equation containing more 
than one degree of derivative, accuracy can be gained for only one 
specified derivative at a time. The MOD 1 approximation has the 
accuracy gain in the first derivative, and the MOD 2 approximation 
has the accuracy gain in the second derivative. 

Examination of Eqs. 3.4, 3.10, and 3.14 reveals that the 

only difference in these three equations is in the third term. 
Y 541 and Yye1? while Eq. 3.4 
does not. This ability to let y vary over the interval j+1 to j-1 is 


Equations 3.10 and 3.14 allow values for 


the reason truncation error is decreased, Jf the variation of y 
happens to be linear over the three successive stations j+l, j, and 
j-l1, Eqs. 3.10 and 3.14 then both reduce to the conventional equation, 
Bq. 3.4. 
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CHAPTER IV 


COMPARISON OF THE SOLUTIONS 


The three finite difference approximations to the differen- 


tial equation, Eq. 3.3, are repeated here for convenience. 








Conventional 
ei aie _c —_ 
Ate \U+1 2@Y, F9y-, } + sor Ue a | 
ame mM (3.4) 
MOD 2 
“ ly peut + Ly g | 
At \U#H i per Z2AtT\Uti U-1 
k o 
+ FE(Ye, TOY + Y, | = (3.10) 
MOD 1 
m 


Jiri 24, TY, |+ sat (Ya Ue! | 


k = 
+ eat oy ty) = a 


The exact solution to the differential equation is of the 





AP 


form 
y=a eG (4.1) 
where A is an arbitrary constant. Substituting Eq. 4.1 into the 
differential equation, Eq. 3.3, leads to 
| mar+ ca + k | 4 et? =O 
If Eq. 4.1 is to be valid for all values of t, then 


raced . =o 


Thus, 
Ee t Vice -—4mk | 
et — 


Zin (4.2) 


20 





If w= k and C = E , then Eq. 4.2 becomes 
Yom 2mw 
a = —wl ei) Vi-¢é 
(4.3) 


C is the nondimensional damping coefficient. A value of 
C less than 1.0 means the system is underdamped, and the motion is 
oscillatory. A value of C greater than 1.0 means the system is 
overdamped, and the motion is non-oscillatory. When is exactly 
equal to 1.0, the system is critically damped. For this analysis, 
the special case of critical damping is not considered. 


The solution assumed for the finite difference equations ae 


J) =A Nod 
(4.4) 
where j is the station number. Solutions for x that satisfy 
Eqs. 3.4, 3.10, and 3.14 are derived in Appendix D. For the first 


order finite difference equation, 


(1-8) 4 jarfi— c?-2 ao? 
. = —Se 4 ee. 
ie a. | 
(4.5) 
where () = Wht. 
A discussion of G) is in order here. The quantity q) is 
the natural frequency of the oscillator in radians per second. 
Therefore, @)/\ft or G) has the units of radians per time interval and 
is a measure of the resolution of the solution. A small value of 
G) would mean that the system is being sampled many times per cycle. 
For example, a value of () equal to 0.5 would mean the system was 
being sampled every 0.5 radians, or over 12 times per cycle. A value of 
G) over 1.0 would mean that the solution has a very poor resolution, 
and the truncation error would be high. 
Equation 4,1 can be written in the form 


J =A @ a (UAt| 


At 
y =A| 2° |? (ane 
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or 


Comparing Eqs. 4.4 and 4.6 reveals that 


Gar ov 


is a proper comparison. However, Eqs. 4.3 and 4.5 show that both 


Cc 


CQ and r are complex numbers for the underdamped system. There- 


fore, 
US Gop ot | Clee 
Cc Cm! C 7 
(4.7) 
where 
ae ae ae: 
a. = w C » 2, =wVI-C 


and for the first order difference scheme 


where 
| io = | = . 2 Pee 
gee | , — NEE 
a= = and — = 
© mee! &) 5 1+l@ 
The solution x can also be given in the form 
(=rel 
where 
=\ / 2 2 = z 
PHENOM + 2 and 4 =ARC SIN = 
A further simplification can be made by replacing f with 
e. Thus, 
K=- In &£ 
and r becomes _ + iO (4.8) 


WS Ge C 
Examining Eq. 4.7 and 4.8 reveals that 


a Aft VS -K 
F 


and 


a, At Vs 6 


are proper comparisons. Let the relative error in both K _ and a 
be defined by 





=~ ea x and V = ae + 
P a. At GO 


Applying the value of Qi and Gd -. to {L and V , leads to 
@ fy - ¢7 C @ 
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The relative error terms, [L and PY , provide an indication 
of the accuracy in the numerical solution. This can be seen by com- 
paring the complete solutions. Substituting Eq. 4.3 into Eq. 4.1 
yields the exact solution 

m= + / To | 
Y=A ol Go | Gade C —_ 
(4.10) 


sincef=jJ/\+. Substituting Eq. 4.8 into Eq. 4,4 yields the solu- 
tions for the finite difference approximations 
| -K +i Oly 
yY=Ac 
(4.11) 
Applying the relative error terms ee and Y given by Eq. 4.9 to 
Eq. 4.11, the finite difference solutions take the form 
[-vew + ipwrT- elo 
Go = Ane 
(4.12) 

Thus, a comparison of Eq. 4.12 with 4.10 shows that if UL and VY 
are both equal to 1.0, Eq. 4.12 is the exact solution. 

Examining Eq. 4.12 reveals that UL and TY have a specific 
role in the error analysis. In Eq. 4.12, -vCwis a measure of the 
rate of attenuation in the finite difference solution. Likewise, 
LU} } - 4is a measure of the frequency of the oscillations in the 
finite difference solution. Therefore, JY will be an indication of 
the error in attenuation, and ML wild be an indication of the error 
in frequency of the numerical solution. 

6 For the overdamped problem o anes not change. However, 
e! 


becomes an attenuating factor and Eq. 4.12 becomes 
= + 2 t 
[-vow thw Py jua 
Y=Ae : C-! 
where PY is given in Eq. 4.9, and JL becomes 


_ aa (4.13) 


Y* OVET 
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The expressions for Ds ae 6 » and 1G for Eqs. 3.4, 3.10, and 
3.14 are derived in Appendix D. Table I lists the values of ) , 

K , and 6 , and Table II lists the values for /@ for all three 
numerical approximations. The subscripts c, 1, and 2 refer to conven- 
tional, MOD 1, and MOD 2 respectively. 

In conclusion, a proper comparison of the solutions to the 
three finite difference approximations with the exact solution has 


been shown to take the form 


V = and YW = <a 
78) ON be 


for the underdamped problem and 


Vege oe 6 


Go and w= Zs 


2 
for the overdamped problem. A computer program Sate the 
necessary comparisons was written in the FORTRAN IV language for the 
IBM 360 digital computer. The results obtained from this program 
follow. 
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TABLE I 


LIST OF THE UNDERDAMPED SOLUTIONS 
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TABLE II 


LIST OF THE OVERDAMPED SOLUTIONS 


/Q. = ARC SINH # 





/@ = ARC SINH +f | 


/ 6, = ARC SINH ?*? 
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CHAPTER V 


RESULTS AND CONCLUSIONS 


Results 
Figure 2 is a plot of the data in Table III. For this plot, 
( is equal to 0.0, and the system is undamped. The undamped solu- 
tion will have no attenuation and therefore, the VY term is meaning- 
less. 


With C = 0.0, the differential equation becomes 


fit) ea = 0 /- 


Recall that the MOD 2 method was derived by maintaining the Hermitian 
form of the second derivative. As a result, the MOD 2 method pro- 
duces very little error for all values of G) . Since the first de- 
rivative is not present in Eq. 5.1, the MOD 1 method can be expected 
to have more error than MOD 2. However, it is found that the error 
in MOD 1 is approximately the same as the error associated with the 
conventional method. | 

The significant gain in accuracy by the employment of MOD 2 
means a larger time interval could be used to approximate the system 
without loss of accuracy. For example, if the natural frequency of 
the system equaled 180 radians per second, a time interval of 0.005 
seconds would describe the system with 0.2 per cent error in 
using MOD 2. To achieve this accuracy with the conventional method, 
a time interval of 0.0015 seconds must be used. This means that the 
time interval can be increased by a factor of three by using the 
MOD 2 approximation for this problem. 

As the damping is increased, the coefficient of the first 
derivative increases, and the first derivative begins to have a sig- 
nificant influence upon the response. The MOD 1 method should become 
more accurate as C increases since it is derived by maintaining the 
Hermitian form of the first derivative, Figures 3 and 4 are plots of 
the data presented in Table iV with C equal to 0.9, Since { is 


less than 1.0, the underdamped comparisons are used. As anticipated, 
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TABLE IIT 


RELATIVE ERROR FOR a 0.0 





wo LHe <3 
{0.99958 [1.00000 _ | 
a 
| 0.98907 | 











| 0.8 41.02879 | 0.97513 [1.00088 _| 
0.9 $1.03726 | 0.96907 f 1.00142 __ 


__ 


q 
figs a, 04720 . 0. 962551 00219 


60 


6 


SNIdGWVG YOS NOILNIOSSY ONMdWVWS sA7? ¢ 3yNdlS 
IWAYSLNI SWIL/QVY VM 
Ke . 3 


G v £ co 











| GOW 


on 





30 


60= ONIdNVG YOS NOILMIOSSY ONIIdWVS $4 4 Bb JAYNOlS 


“IWAYSLNI SWIL/QV8 VM 





ol 


€ | 


v | 


S| 


9] 


31 


TABLE IV 


RELATIVE ERROR FOR C = 0,9 
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the MOD 1 approximation is seen to be the most accurate, For example, 
in Figure 3 at a resolution of 0.8 the MOD 1 approximation has 3 per 
cent error in 6 . For this same resolution, the conventional approxi- 
mation has 42 per cent relative in 6 ; 

Figures 5 and 6 are plots of the relative error versus C 
for a constant resolution of 0.5 as listed in Table V. Figure 5 
shows that one of the two Hermitian approximations always has a better 
value of ee than the conventional approximation. For low values 
of C the MOD 2 approximation is best, and for high values of iG the 
MOD 1 approximation is best. This means that throughout the range 
of C , greater accuracy in the frequency can be obtained by using one 
of the Hermitian approximations. 

The error in attenuation, ) , is shown in Figure 6. Note that 
for the lightly damped system, where the attenuation is small, the 
conventional approximation has the least error. However, this does not 
imply that the total solution given by Eq. 4.12 is most accurate when 
the conventional approximation is employed since Figure 5 shows that 
the conventional approximation has the largest error in LM for the 
lightly damped system. For values of C between 0.45 and 0.65, the 
MOD 2 approximation is best, and for values of C greater than 0.65, 
the MOD 1 approximation is best. Thus, Figure 6 shows that the Her- 
mitian approximations have less error in attenuation for values of ( 
greater than 0.45. 

Analysis was not made for the critically damped oscillator. 
However, as C approaches 1.0 from either side, Figure 5 shows that 

ML becomes very large. In order to understand what happens to the 
numerical solution for .9< C< 1.1, it is necessary to make other 
comparisons. For example, the exact solution, Eq. 4.10, and the 
solution for the finite difference approximations, Eq. 4.12, can both 


be given in the following form 


y= 4 2[-§ + ij ViI- CO) De ol 
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TABLE V 
RELATIVE ERROR FOR W = 0.5 






—s ———— oe 


Pa : 
ae 1.00013 | 


|i. 00074 


ewer (1.01531 | .9666| 0. 99598 | 9873 |. 00551 § 
1.0140 §1.01849 | .9725 j 1.00070 ° 


— san, BA i i ee Og St 











1.2101 | 1.28394 | 1.2517 | 1.42483 












p15 1 1.2970 { 1.40661 
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for CE 1.0. Thus a plot of,// — & andi) i= ee versus C shows 


the error in the oscillatory part of the solution for the underdamped 
system. For the overdamped system, C 7 i.0, Eqs. 4.10 and 4.12 


is, id -6 ¢ VT] wate 


and [-vOCtHNVG-] lwaty 


Y=Ac 


Thus, for this range of C a plot of (7, ane ene C 


will be the equivalent comparison for the overdamped system. 
Examination of the equations in Table I reveals that the 
numerical solution will be critically damped for the values of C slightly 
less than 1.0. Therefore, to show the critical damping point for each 
finite difference approximation it is necessary to use the comparisons 
for the overdamped system for values of .95€CK1.0, Equation 4.13 
defines {1 for the overdamped system as 
Lo ee 
WW 2 
mm 
/ 0 


tN te ee 


therefore, 


(5.2) 
Using Eq. 5.2, it is possible to evaluate the overdamped comparisons 
for values of e- 1.0. 

Figure 7 is a plot ofl] - for LV CE [fox values of G 
between 0.9 and 1.1 with W = 0.5. The MOD 1 approximation has its 
critical damping point closest to 1.0, and the conventional approxima- 
tion has its critical damping point furthest from 1.0. As W) is de- 
creased and the system is sampled more often, the critical damping 
points of all the numerical solutions will approach 1.0. However, the 
Hermitian approximations will always be critically damped closer to 
1.0 than the conventional approximation. The fact that' each numerical 
solution has its own critical damping point that is less than 1.0 
explains the large values of [ML in this region. 

For the overdamped system a problem arises when CG is greater 
than 1.0. Examination of the equations in Tables I and II reveals an 
unstable zone for these finite difference equations at this point. 


Although ( G@) = 1.0 is essentially the limit for numerical stability, 
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the accuracy of the difference approximations becomes very poor 
before the unstable zone is reached. As the stability limit is 
approached, both LL and YY approach infinity. Figure 8 shows the 
stability limits on all three finite difference approximations. This 
plot was made using the equations in Tables I and II. 


Conclusions 

From the results presented here, the conclusion can be made 
that the Hermitian equations can be used to obtain more accurate 
numerical solutions for the same size increment than can be obtained 
by first order central finite difference equations when applied to 
differential equations containing derivatives of more than one degree, 
Further, this increase in accuracy is not accompanied either by an 
increase in complexity of the resulting matrix when applied to boundary 
value problems or by a significant increase in computation time when 
applied to initial value problems. Consequently, use of the Hermitian 
equations would be most advantageous when computer core storage space 


is limited and accurate results are required. 
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APPENDIX A 


THE DERIVATION OF THE CONVENTIONAL EQUATIONS 
FOR THE FIRST AND SECOND DERIVATIVE 


The first order central finite difference equations can be de- 
rived from the Taylor series. The stations are lableled station 
jtl, station j, and station j-l. The interval between stations j+1 
and j is called fA\ . 

The Taylor series for Y j41 is 


Fi 2 aa 
= 4 1A i he, 7 pay Iv 
Fe J, t AY Os) eid) eo + 
o 
ryt 
/ (A.1) 
For Yeo. the Taylor series becomes 
= _ / i YW A? Stl 
je) wd} Be 2) ey 
2a aye se 
+, Y- ery * (A.2) 
Subtracting Eq. A.1 from Eq. A.2 gives 
‘. D2”, Bi ,... 
= = -~- => aa r 
Yer J 72754 * TY * coy 
(A.3) 
Solving for y'; leads to 
ee _ 7 A2, IS Kz _* 
Y= TK (Ya id 6g ney” he 


This is the first order central finite difference equation for the 
first derivative. 


Adding Eqs. A.1 and A.2 gives 





tt / At IZ 
i i ee | Yi)~ BY, 
LNG v7 @ e @ 
~ 3604, 7 (A.5) 


Equation a.5 is the first order central finite difference equation 


for the second derivative. 
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The truncation error in each equation is usually approximated by 
the first term in the Taylor series that is not used in the difference 
approximation. For example, the error in Eq. A.4 is 

Xe Wt 
“6 wv 
and the error in Eq. A.5 is 


TON Zz Iv 


VD 4) 
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APPENDIX B 


THE DERIVATION OF THE HERMITIAN EQUATIONS 
FOR THE FIRST AND SECOND DERIVATIVES 


Taking the derivative of Eq. A.1l in Appendix A gives 
2 $ 
Wj Ly MW A 
— + —_a_vE 
At. .£ AF a 
————e - — (a Jd 
tre Yt 20% 


/ it 


Y= 


Likewise Eq. A.2 becomes 
/ / tw f2 mW AF 
Yer = YO OY THy- ey 
At z A? 
‘IE 4 7704 
Adding Eqs. B.1 and B.2 gives 


ft 74° 


/ / a / as “U1 (aS Be ae, 
Fi 4/ . 4-1 =AY Fr 8 J) " 72 f i 
Dividing Eq. B.3 by 6 and solving for Ne y''' leads to 
6 
A2 Ml _ / / / / _ Os Z -_ 
SY = ee 24" Yul ae? 


From Eq. A.4 of Appendix A, 





i. | oa £.5,7 
4 7A We 4y-,]- SY) - jo 
Thus, substituting Eq. B.4 into Eq. A.4 gives 
ee = _ ef / / 
Y= FAlYya y.| oY, | me 
a ae | 
+ 7904, ae 


44 


(B.1) 


(B.2) 


(B.3) 


(B.4) 


(A.4) 


(B.5) 








Solving for Y iad - Y 5 oh gives 


= mye / / / 
Yor ~ Yor = BY AY Ye 
ae ce 

79 VU (B.6) 


Equation B.6 is the Hermitian equation for the first derivative. 


The second derivative of equation A.1 is 
W HW um we ew Lf £ 
— ._: —— 
7 them th” > Uemilaetiaieiaas ae 
ip 5 saT 
a Vib 
ba, O* 
24 UY ~ 120 


o 6 @ 








(B.7) 
Likewise Eq. A.2 becomes 
: Ae Ae 
Yi “YO AW" TH EY 
+ 5 we 
+ 4Liyt Dy 4472 @ 
AFUY (207 (B.8) 
Adding Eqs. B.7 and B.8 gives 
14 yf // 2 Ww A? WL 
+ = : + &y +-°° 
a) yn oy “ny 2 
(B.9) 
Dividing Eq. B.9 by 12 and solving for A? y ™ yieids 
eZ 
A2E_ yu ” v\ At 
i272) ~ 72 Yiu 2Y t Ziljey terre 
(B.10) 
From Eq. A.5 of Appendix A, 
ji Mi ots v4 
d= AY 24+ YY, 
STE 
“See 7s (A.5) 
Substituting Eq. B.10 into Eq. A.5 gives 
= 1 = IG] + Yy )- ial; 7 + y” 
d) SRS A a ia are 0 iy, Ua! 
¥ a Latins. 
Z#O 4 
(B.11) 
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Then, solving for (Yuua - ayy + Y yy) gives 
> IN M4 Ne WT] // 
a - $ i 
Ye 2 WP 24 * Y, 
6 





Equation B.12 is the Hermitian equation for the second derivative. 


A comparison of the conventional equations in Appendix A with 
the Hermitian equations derived here reveals that the Hermitian 
equations have less truncation error. This is to be expected since 


the error term in each first order equation was used to derive the 


respective Hermitian equation. 
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APPENDIX C 


DISCUSSION OF THE TRUNCATION ERROR 
IN THE MODIFIED HERMITIAN EQUATION 


In the MOD 2 method of analysis, the Hermitian equation for 
the first derivative is modified. Equation B.6 is the Hermitian 


equation for the first derivative 














oy =-Aly’ <9 eae 
4 Ye =F (Yur *Y 6 Yd” FoY 
(B.6) 
This equation can be written 
Cc — / / / 
max | Mere Y) _)= — Wn t TY F Yr, | 
_ A? yz 
360 U (ca) 
The conventional equation for y' j can be written : 
Cc c AZ 17 
_—_ ~ —— a ie 
ae Yi TY te 
(C.2) 
Addition of equations C.1 and C.2 gives 
mis i aes ie 
2A IY). y_)= (2 We as a, 
by ee 
ie J 369~) 8 


Equation C.3 is the modified Hermitian equation for the first 
derivative that was used in the MOD 2 approximation. The truncation 


error in this equation is 
A+ stl 
72 


This truncation error is one-half the truncation error in the first 





order equation for the first derivative. Therefore; in the MOD 2 
approximation to the differential equation, there is very little gain 


in accuracy associated with the first derivative. 
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Likewise, for the MOD 1 approximation, the Hermitian equation 


for the second derivative must be modified. The modified equation is 


ym . _™m 4¢ 1 LA 
Bly 2% * Yu E Yt ty Ye 


(C.4) 
The truncation error, 
A*-* ¢£ 
72% 
is equal to the truncation error of the first order difference 
equation for the second derivative. Therefore, in the MOD 1 approxi- 
mation, there is no gain in accuracy associated with the second 


derivative. 
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APPENDIX D 


SOLUTIONS FOR THE 
FINITE DIFFERENCE APPROXIMATIONS 


The conventional approximation will be used as an example for 
this appendix. The conventional finite difference expression for 


the differential equation is 


———- 





m i, 
APY 7% * 4.) + Ur Y-1 


+k = O 
Y 0.1) 


When the assumed solution 


J 
yY =A dX 
J 
is substituted into Eq. D.1l, the result is 


wy of 
sje k [AX 








aw [Az 2- Xe] + aoe 
ae 


For an arbitrary solution, 


4 +k). =O 
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(A eee 
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By using the quadratic equation, 








- pt Aap say 
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or [2 m “Kays \ [lam -k (At; 4m? = car) 


i 2| m ~ £Ar| 


When the substitutions 


mg ) C 2 Muy) and. W=W Oe 


are made, Ne becomes 


= = 2 Ee 
5 i -$a7] + @ C =ltee 


A= 


Cc 


| + CB 
To make } ~ complex, take > [-1 outside the radical. Then 


i-$a]eiw Vi — F - Fa 
ee 
, + C@ ial 


The solution IN P is now of the form 
r = x+/2 
Cc 


To get r a into the form 


Y= relG 
G. 
the following relationships are used. 
— 2h me Q@.= ARC SIN 
-— x + Ze and co. (ps 
Therefore, 
2 : Vy 
Eos ers ee te =o 
4 _+{l ta']+@ | - > er 1] 
1+ C¢@ 
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or pox uy ee % 


if, K =-— Io r, then 


K= Irn tn (alae, G) 
l + (w (D.3) 


Similarly, 


0. =ARC SIN = 





(D.4) 
For the overdamped system, i OH. must be found, From Eq. D.4, 





Therefore, 


/@ = ARC SINH? 





(D.5) 
The equations for » kG <5 6 » and if for MOD 1 and MOD 2 


are derived in a similar manner. 
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